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ABSTRACT 



We present the sub-mm/mm SED for a sample of eight young circumstellar disks in the outer regions of the Orion 
Nebula Cluster. New observations were carried out at 2.9 mm with the CARMA array and for one disk, 216-0939, at 
3.3 and 6.8 mm with ATCA. By combining these new millimeter data with literature measurements at sub-millimeter 
wavelengths we investigate grain growth and measure the dust mass in protoplanetary disks in the Orion Nebula 
Cluster. These data provide evidence for dust grain growth to at least millimeter-sizes for the first time in a high-mass 
star forming region. The obtained range in sub-mm/mm spectral index, namely 1.5-3.2, indicates that for disks in the 
outskirts of the Orion Nebula Cluster (projected distance from the cluster center between about 0.4 pc and 1.5 pc) 
grain growth to mm sizes occurs in the same manner as disks in regions where only low-mass stars form. Finally, in 
our sample three disks are more massive than about 0.05 M©, confirming that massive disks are present in the outer 
regions of the Orion Nebula. 

Key words, stars: planetary systems: protoplanetary disks — stars: planetary systems: formation — stars: formation 



1. Introduction 

Disks around young Pre-Main Sequence (PMS) stars are 
thought to be the cradles of planets, and for this reason 
they are often called "proto-planetary disks" . Evidence for 
the very early stages of planet formation in protoplanetary 
disks has been obtained by a variety of different techniques, 
ranging from the optical to the millimeter-wave spectral 
domain. Compared to observations at shorter wavelengths 
, the inspection of the sub-mm and mm emission of disks 
. provides two key advantages for the investigation of the 

■ physical processes related to planet formation. First, except 
, for the very inner regions of the disk, the dust optical depth 

■ at these long wavelengths is sufficiently low that the disk 
I midplane can be probed. Second, since the maximum grain 

size to which a particular technique is sensitive to is of the 
order of the observational wavelength, the (sub-)millimeter 
spectral regime offers the chance to investigate grain growth 
up to mm/cm sizes (see Rodmann et al. I2006p . 

So far nearly all the information we have on the dust 
grain properties in protoplanetary disks has been derived 
from systems in low- mass Star Forming Regions (SFRs). 
Little information on the dust grains in young disks is 
present for high-mass SFRs. The analysis of dust in high- 
mass SFRs is particularly compelling since about 75% 
of the observed Young Stellar Objects (YSOs) formed 
in a clustered environment, where high-mass stars form 
(Bressert et aL l^nTU]) . 



Some evidence for dust grain growth in protoplanetary 
disks in the Orion Nebula Cluster (ONC) has been obtained 
at optical and infrared wavelengths. Throop et el. (|200ip 
and Shuping et al. (|2003l) observed in optical and NIR the 
largest ONC silhouette disk (114-426). From the exinction 
curve of the background nebular emission they found that 
the NIR opacity of the disk is dominated by grains as large 
as about 10 fim, larger than the typical sub-micron sized 
grains found in the ISM (Mathis et al. I1977p . Using low- 
resolution spectroscopy in the mid-IR, Shuping et al. (2006|l 
found evidence for grain growth to a few fj,m in the surface 
layers of eight disks in the Trapezium Cluster. 

During the last 15 years several ONC disks have 
been surveyed by interferometric observations at mm- 
wavelengths (e.g. William s et al. ?IUU5[ Eisner et al. [MM 
Mann & Williams I2009b[) . However, the targeted disks were 
quite close to the main ionizing source of the cluster, the 
06-type Trapezium star 8^ Ori C (typical projected dis- 
tances lower than 0.2 pc by assuming a distance to the ONC 
of 40 pc, M enten et al. lMITl Sandstrom et aL lMITl Kraus 
et al. I2009p . As a consequence, the high levels of free- free 
emission by ionized gas swamped the dust emission from 
the disks, making estimates of grain growth almost impos- 
sible to derive from multi-wavelength dust measurements 
in the millimeter spectral window. 

In this work we present new mm-wave interferometric 
observations for a sample of eight disk systems in the out- 
skirts of the ONC. Because of their location relatively far 
from 9^ Ori C (projected distances between 0.4 and 1.5 pc. 
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Fig. 1. Location of the eight disks in our sample (within 
white circles) overlaid on a JCMT-SCUBA 850 /xm image 
of the ONC (Johnstone & Bally [HSni)- The black and white 
star near the bottom of the image marks the position of 
the 9^ Ori C star. White crosses represent the positions of 
proplyds identified by HST observations. 

see Figure [IJ , information on grain growth can be derived 
from mm-wave observations. This allows us to probe the 
dust properties in these ONC disks over the entire sub- 
mm/mni SED. 

2. Observations 

We describe here new mm-wave observations of eight young 
circumstellar disks in the ONC. The disks were selected for 
being detected within the SMA survey of protoplanetary 
disks in Orion (Mann & Williams 2010, in press) with a 
0.88 mm-flux density larger than about 15 mJy, and by be- 
ing located in the outer regions of the ONC (i.e. projected 
distance from the Trapezium Cluster larger than 3'). This 
latter selection criterion was chosen to avoid contamina- 
tion at mm-wavelengths from the non-uniform background 
emission and free-free emission from ionized gas typically 
observed in the inner ONC. The disks in our sample are 
listed in Tabled] 

The dust thermal emission toward our sample of eight 
young disk systems in the ONC was observed with the 
Combined Array for Research in Millimeter Astronomy 
(CARMA) between 2010 Mar 26 and Apr 01. The array 
was in C configuration which provides baselines between 
30 and 350 m. Observations were carried out at a central 
frequency of 102.5 GHz (2.92 mm), with a total bandwidth 
of 4 GHz. 

The raw visibilities for each night were calibrated and 
edited using the MIRIAD software package. Amplitude and 
phase calibration were performed through observations of 
the QSO J0607-085. Passband calibration was obtained by 
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Fig. 2. Flux density at 2.9 mm versus observing day for the 
QSO J0609-085 used as the phase and amplitude calibra- 
tor during the CARMA observations. Red circles represent 
fluxes obtained using Mars as the flux calibrator, whereas 
blue circles have been obtained using Uranus. The dotted 
line represents a linear flt of the data, and the two contin- 
uous lines refer to the ±15%-levels from the linear fit. 



observing the QSO 0423-013. Mars and Uranus were used 
to set the absolute fiux scale. The uncertainty on CARMA 
flux calibration is typically estimated to be ^ 15% and is 
due to uncertainties in the planetary models and in the cor- 
rection for atmospheric effects and hardware instabilities. 
Figure [2] shows the measured 2.9 mm-fiux density of the 
QSO J0607-085 during the period covered by the CARMA 
observations. Both in the case that the QSO flux was slowly 
fading during the week of observations (with a light curve 
shown by the dotted line in the plot, which represents a 
linear fit of the data), and in the case that the QSO had a 
constant 2.9 mm fiux density (equal to the weighted mean 
of the measured fluxes) , the repeatability of our flux mea- 
surements is within the adopted 15%-level. 15% is therefore 
a conservative estimate since part of the scatter may be due 
to intrinsic variability of the QSO on 1-day tiniescale. 

Maps of the dust continuum emission were obtained by 
standard Fourier inversion adopting natural weighting, and 
photometry was obtained in the image plane. The resulting 
FWHM of the synthesized beam is about 2". 

One of our disks, 216-0939, was also observed with the 
Australia Telescope Compact Array (ATCA) and the new 
CABB digital filter bank, which provides a total contin- 
uum bandwidth of 2 GHz. Observations were carried out 
at a central frequency of 91.000 GHz (3.294 mm) on 2009 
Oct 13, and of 44.000 GHz (6.813 mm) on 2009 Oct 14. The 
ATCA array was in the hybrid H168 configuration, provid- 
ing an angular resolution of about 3" at 3.3 mm and 6.3" 
at 6.8 mm. 

The gain was calibrated with frequent observations of 
0529-1-075. The passband was calibrated using 1921-293, 
and the absolute fiux scale was determined through obser- 
vations of Uranus. The uncertainty on the calibrated flux 
is about 30% at 3.3 mm and 20% at 6.8 mm. The MIRIAD 
package was used for visibilities calibration, Fourier inver- 
sion, deconvolution and imaging. 
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Table 1. Sub-Millimeter and Millimeter Flux Densities and Disk Sizes 
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HST Size 
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(arcsec) 


(arcsec) 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


121-1925 


15.0±1.5 


2.4±0.4 






1.5±0.4 


Unresolved 


0.8 


0.003 


0.1 


132-1832 


16.5±1.7 


<2.7 






>1.5 


1.1±0.4 


1.3 






136-1955 


77.6±1.2 


1.8±0.4 






3.2±0.5 


0.5±0.1 




0.02 


0.006 


141-1952 


30.6±1.2 


1.9±0.4 






2.3±0.5 


0.4±0.1 


0.7 


0.02 


0.003 


181-825 


54.8±1.0 


4.4±0.5 






2.1±0.4 


Unresolved 


0.7 


0.01 


0.05 


216-0939 


91.9±1.3 


3.5±0.6 


2.5±0.3 


0.34±0.05 


2.7±0.2 


1.6±0.1 


2.6 


0.04 


0.02 


253-1536^ 


134.2±1.0 


7.9±0.7 






2.6±0.4 


1.5±0.1 


1.5 


0.04 


0.05 


280-1720 


21.8±0.7 


<1.8 






>2.1 


Unresolved 


0.8 







1) Object name. 2) 0.88 mm-flux density (from Mann & Williams 2010, in press). 3) 2.9 mm-flux density. 4) 3.3 mm-flux density. 
5) 6.8 mm-flux density. 6) Sub-mm/mm spectral index derived by fitting with a power-law the data shown in the table. The 
errors account for the uncertainty on the flux scale at each waveband. 7) FWHM along the disk major axis from a gaussian fit in 
the SMA images. 8) Disk diameter from HST images (Smith et al. .2005t Ricci et al. ;2008). 9) 2.9 mm-dust opacity derived by the 
SED-fitting procedure (and the adopted dust model) assuming q = 2.5. 10) Disk mass derived by the SED-fitting procedure 
assuming q = 2.5 (see Sect. 13.3")) . a) The SMA observations of 253-1536 show a weaker secondary component 
(-Fo.ssmm = 37.6 ± 1.0 mJy). The spectral index of this system has been calculated by using for Fo.ssmm the sum of the SMA 
fluxes separately measured for the two circumstellar disks. For this reason the derived spectral index and disk mass are 
descriptive of the whole system rather than of the two disks separately. 



3. Results 

Table [1] reports the sub-mm/mm flux densities obtained 
for our sample of eight young circumstellar disks in the 
outskirts of the ONC. Column (6) lists the sub-mm/mm 
spectral index a {F^, ^ j/") as derived from the results of 
our observations. Whereas all the eight disks in our sam- 
ple have been detected with the SMA at 0.88 mm (Mann &; 
Williams 2010, in press), we detected six of them at 2.9 mm 
with CARMA. For the two undetected disks, 132-1832 and 
280-1720, the upper-limits in the 2.9 mm-flux density trans- 
late into lower-limits for the sub-mm/mm spectral index. 
The range in sub-mm/mm spectral index spanned by our 
sample, i.e. ~ 1.5-3.2, is very close to what observed for 
class II YSOs in low- mass SFRs (see Figure [3]). This in- 
dicates that the properties of dust grains for disks in the 
outskirts of the ONC are very similar to those found in 
disks in low-mass SFRs. Furthermore we did not find any 
clear trend between spectral index a, HST size, disk mass, 
or distance from 6^ Ori C. 



3.1. Grain Growth 

Information on the grain-size distribution of dust in the 
outer regions of protoplanetary disks can be obtained by 
analyzing the disk SED at sub-mm/mm wavelengths (e.g. 
Draine [2006p . More specifically, for a disk whose sub/mm- 
mm flux is entirely due to optically thin dust emission in 
the Rayleigh- Jeans regime, a very simple relation between 
the measurable spectral index a of the sub-mm/mm SED 
and the spectral index /? of the dust opacitjQ holds, namely 
/3 = a - 2. 

Values of /? lower than the value found for the unpro- 
cessed ISM dust population (/3ism ^ 1.7-2) can be obtained 
only with dust population where grains as large as at least 
~ 1 mm are present (see e.g. Natta et al. I2007p . Deviations 
from the aforementioned relation between a and /3 can be 
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Fig. 3. Spectral index between 0.88 and 2.9 mm plotted 
against the flux at 0.88 mm. Black stars and lower-limits are 
for the ONC disks, whereas grey circles and squares are for 
class II disk s in Tau rus (Ricci et al. I2010ap and Ophiuchus 
(Ricci et al. I2010b[l , respectively. The fluxes of the Taurus 
and Ophiuchus disks have been scaled to the ONC-distance 
adopted in this work (400 pc) to allow a comparison with 
the ONC disks. In the lower-right corner of the plot, the 
typical errorbars of the plotted quantities are shown. 



accounted for by using physical disk models to fit the disk 
SED together with the knowledge of the spatial extension 
of the disk, which is needed to estimate the impact of the 
inner optically thick regions on the long- wave spectrum (see 
Testi et al. I200ip . By doing this for 38 class II disks around 
low-mass PMS stars in the Taurus- Auriga and Ophiuchus 
SFRs, Ricci et al. (|2010a[ 2010b) found a-P = 1.8±0.1. AU 
disks in our sample have outer radii larger than ~ 100 AU 
(Table [1] col. 7 and 8), all the model fits have only a mini- 
mal contribution from optically thick emission (see discus- 
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sion in Ricci et al. l2010a| . Applying the a-/3 relation to the 
Orion disks, we obtain values for the dust opacity spectral 
index fi ranging between about for 121-1925 and 1.4 for 
136-1955. Apart for the case of 136-1955 whose Icr-interval 
of P is consistent with /3ism, for all the other five disks 
detected with CARMA the constrained /3- values are signif- 
icantly lower than /3ism, and thus we found firm evidence 
of grain growth to '~ mm-sized pebbles for the first time in 
the ONC, and, more generally, in the cluster environment 
of a high-mass SFR. 

Note that, to derive the results described so far, we as- 
sumed that all the measured sub-mm/mm fiux comes from 
dust emission only. At mm- wavelengths some contribution 
to the measured flux from ionized gas has been found in the 
literature. This contribution can be constrained from ob- 
servations at longer cm- wavelengths, where free-free emis- 
sion from ionized gas dominates over the dust emission. For 
disks in low-mass SFRs the presence of ionized gas is prob- 
ably induced by the ionizating flux from the central star, 
and it has been found to contribute typically for 25% or 
less at 7 mm and a few percent at most at 3 mm (Wilner 
et al. [MI^ Rodmann et al. [200^ Lommen et al. 120(1]) ■ 
For disks in the inner ONC at angular separations lower 
than about 1 arcmin from the main ionizing source of the 
entire region, 9^ Ori C (projected separations lower than 
about 0.1 pc), free- free emission has been typically found 
to dominate the emission at wavelengths close to 3 mm 
and longer (Williams et al. 120051 Eisner et al. 120081 Mann 
& Williams I2009bp . In order to get accurate estimates for 
dust grain growth from multi-wavelengths observations in 
the millimeter, it is important to consider disks which are 
further away from 6*^ Ori C, for which the free-free emis- 
sion does not swamp the mm-wave emission from dust. Our 
sample is made of sources which are at projected distances 
between about 0.4 and 1.5 pc. Since our disks are much fur- 
ther away in terms of projected distance to 9^ Ori C than 
the disks observed in the literature, in the rest of the paper 
we will assume that no contribution from free- free emission 
is present for the disks in our sample. This is supported by 
the fact that the derived spectral slope a for our sources do 
not correlate with the projected distance to 9^ Ori C, as one 
may expect if significant contamination from free- free emis- 
sion was present in our data. Furthermore, in the case of 
216-0939 for which new data at 3 mm and 7 mm have been 
obtained, the sub-mm/mm SED does not show any change 
of the spectral slope (Figure S]) , as it would be expected if 
free-free emission was significant at these wavelenghts (see 
e.g. Williams et al. 12005^ . However, future observations at 
longer cm-wavelengths are needed to really constrain the 
impact of free-free emission and test this assumption. 

3.2. The case of 216-0939 

Figure 2] shows the fit of the sub-mm/mm SED of 216- 
0939, which is the best characterized disk in our sample, 
with two-layer (surface+midplane) models of flared disks 
heated by the radiation of the central star (Chiang &; 
Goldreich lTWl DuUemond et al. 2001). For the steUar pa- 
rameters we considered a PMS star with a K5-spectral type 
(from Hillenbrand 119971) , which we converted into an effec- 
tive temperature of 4350 K using the temperature scale of 
Schmidt-Kaler (1982,) . We then derived a stellar mass of 
about 1.2 Mq and a stellar luminosity of 2.5 Lq by plac- 
ing the PMS star onto a 1 Myr-isochrone from the Palla & 



Stabler (|1999p PMS evolutionary models. For the disk we 
adopted a truncated power-law mass surface density with 
an exponent of 1 (as found by Andrews & Williams 120071 
from disks in Taurus and Ophiuchus) and an outer radius 
of about 290 AU from the size of the mm-emission in the 
SMA map (Table [1]). A disk inclination of 75° has been 
taken from an analysis of optical HST images (Smith et 
al. 2005). From the SED- fitting procedure (Fig. ij) we de- 
rived (3 = 1.0 ± 0.3, which is indicative of grain growth in 
the 216-0939 disk. 

By using the same dust model discussed in Ricci et 
al. ([2010a[ 2010b), i.e. spherical porous grains made of sili- 
cates, carbonaceous materials and water ice (abundance of 
each species from Semenov et al. I2003|) . we could get infor- 
mation on the grain-size distribution, here assumed to be 
a truncated power-law with index q and maximum grain 
size Oniax- A good fit of the SED could be obtained only 
for q < 3.5. Grain-size distributions with (/-values of 3 and 
2.5 fit equally well the sub-mm/mm SED but with different 
flmax-values, namely about 2.5 and 2 mm, respectively, but 
here it is important to note that the precise estimate for 
the maximum grain size Omax can strongly depend on the 
adopted dust model. 

3.3. Derived disk masses 

Other than providing information on the grain-size distri- 
bution, the SED-fitting procedure described in the last sec- 
tion returns estimates for the dust mass in the disk, and, 
by assuming the ISM- value of about 100 for the gas-to-dust 
mass ratio, for the total disk mass. Column (9) in Table [T] 
reports the derived disk masses for the six disks detected at 
millimeter wavelengths, for which information on the dust 
properties could be inferred. For these estimates we consid- 
ered a g— value of 2.5 because for the two disks showing the 
lowest mm-spectral indeces, 121-1925 and 181-825, models 
with g <; 3 cannot reproduce the data. This shows that the 
grain size distribution in these disks is shallower than found 
for the ISM (^ism ~ 3.5, Mathis et ah IW?!) . This resuh, 
which, as described in Ricci et al. (|2010ap does not depend 
on the adopted dust model, is another indicator for dust 
processing in these disks. 

The range in disk mass covered by our sample is 
0.003—0.1 Mq. This shows that massive disks are present 
in the outer regions of the ONC, confirming recent results 
by Mann & WiUiams (12009a', 2010, in press). Note that 
121-1925 is the faintest disk at 0.88 mm but is nevertheless 
detected at 0.88 mm. We thus derive a very low value for 
a. In our analysis, this implies a very low value for f3, and 
a very significant inferred grain growth, suggesting a very 
low value for the dust opacity coefficient at 2.9 mm (see 
col. 9 in Table [J), and consequently a large disk mass. 

4. Summary 

We present new millimeter interferometric observations of 
eight disks in the outskirts of the ONC. Together with sub- 
mm data in the literature, these new data allow us to study 
dust properties, like grain growth and disk mass, from the 
entire sub-mm/mm SED. Except for one disk, 136-1955, 
whose emission is consistent within the errorbars with a 
ISM-like dust population, we found evidence, for the first 
time in cluster environment in a high-mass SFR, of grain 
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Fig. 4. Fit of the sub-mm/mm SED of the 216-0939 Orion 
disk with a two-layer passive flared disk modeL The param- 
eters of the model shown here are described in Sect. 13.21 The 
uncertainties in the data points account for the uncertainty 
on the flux-scale. The color image of the 216-0939 disk is 
from HST observations by Smith et al. (|2005p . 




growth to at least mm-sized pebbles for 5 out of the 6 de- 
tected disks at mm-wavelengths. As for the disk masses, 
estimated by using a dust model which accounts for grain 
growth in the calculation of the dust opacities, our results 
conflrm that massive disks (M > 0.05 Mq) can be found in 
the outskirts of the ONC. 
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